



'/ *,~ r I 'l' t 



J!W’)!* 









P.& 



RAPID APPROXIMATION OF ' 
SERVOMECHANISM FREQUENCY RESPONSE 

WILLIAM F. SALLADA . .I 

and 

VAUSHN E. WILSON, JR. 






K?" J*r^ - '- -- '-fc 
r>.w^ : 



NPS ARCHIVE 
1962 

SALLADA, W. 












library 



DUDLEY KNOX uiBRARY ^ 

NAVAL POSTGRADUATE oCHOO* 
MONTEREY CA 33943-5101 



POSTGRADUATE SCHOOL 
''.‘‘rnRrT!A 



I 



RAPID APPROXIMATION OF SERVOMECHANISM 



FREQUENCY RESPONSE 



■k-fckiri' 



William F, Sallada 



and 



Vaughn E. Wilson, Jr. 



RAPID APPROXIMATION OF SERVOMECHANISM 



FREQUENCY RESPONSE 



by 

William F. Sallada 

Lieutenant Commander, United States Navy 

and 

Vaughn E. Wilson, Jr. 
Lieutenant, United States Navy 



Submitted in partial fulfillment of 
the requirements for the degree of 

MASTER OF SCIENCE 
IN 

ELECTRICAL ENGINEERING 



United States Naval Postgraduate School 
Monterey, California 



1962 



nPc) A rcV\v^ 

j kjO 



LIBRARY 

U.S. NAVAL POSTGRADLfATE SCHOOL 
' MONTEREY, CALIFORNIA 

RAPID APPROXIMATION OF SERVOMECHANISM 
FREQUENCY RESPONSE 



by 

William F. Sallada 
and 

Vaughn E. Wilson, Jr. 



This work is accepted as fulfilling 
the thesis requirements for the degree of 

MASTER OF SCIENCE 
IN 

ELECTRICAL ENGINEERING 
from the 

United States Naval Postgraduate School 



ABSTRACT 



The purpose of this thesis was to develop a technique by which the 
frequency response of a servomechanism might be rapidly evaluated by 
graphical methods. The approach consisted of developing several frequency 
response parameters related to the location of the poles and zeros of the 
closed loop transfer function of a servomechanism. The systems considered 
were the pure second order (one pair of complex conjugate poles) , second 
order with one zero (one pair of complex conjugate poles and one real 
zero), and a pure third order (one pair of complex conjugate poles and 
one real pole). For the latter two systems, the complex poles and the 
frequency were redefined as referred values relative to the real pole 
or zero. 

Loci of constant frequency response parameters were mapped onto a 
complex plane; or loci of constant imaginary referred pole locations were 
mapped onto an auxiliary complex plane, depending upon the specific 
system considered. Using the points obtained from the appropriate 
relationship, the frequency response, M and N contours, could be rapidly 
and easily sketched. The M and N contours could then be used as a basis 
for rapid system evaluation in the frequency domain. 

The authors wish to thank Professor George J. Thaler for bringing 
to their attention the thesis by Major Walter B. Patton, USMC , and 
LT William B. Abbott, III, USN, without which this thesis topic would 
probably not have been considered /!/. 

The authors are deeply indebted to Assistant Professor Robert D. Strum, 
who acted as their thesis advisor, for the friendly support, encour- 
agement, and advice given so selflessly. 
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1. INTRODUCTION 



Each year the Importance and utilization of servomechanisms and 
automatic control systems increase. It naturally follows that more 
papers are submitted pertaining to the design, analysis, synthesis, 
and compensation of servomechanisms and automatic control systems. 

Many of these papers and existing texts discuss details and laborious 
calculations which are certainly a prerequisite to analyze or synthesize 
the operation of a servomechanism. However, in many cases, a rapid 
approximation of system operation would be beneficial, and possibly 
only an approximation would be required upon which an Intelligent 
estimate could be made on the performance of the system under study. 

There are only a few papers known to the authors in which the use of 
approximations are presented from which an estimate of system performance 
can be made, e.g., Floyd's method /2/. One such paper by Patton and 
Abbott /!/ presents a method for a rapid approximation of servomechanism 
transient response. 

A corollary should exist to the Pat ton -Abbott development in which 
a rapid approximation of the frequency response of a servo could be 
determined. The use of frequency response analysis is still an Important 
criterion in measuring servo performance, although the transient response 
approach has become very popular and widely used in recent years, almost 
to the point of neglecting frequency response. There are many ways of 
determining frequency response from transient response, but lengthy 
calculations are necessary. 

It is the purpose of this paper to present a graphical means for 
determining a rapid approximation of servomechanism frequency response. 
This approach eliminates the use of the Bode diagram and Nichols' 
chart. The M and N contours will be obtained directly from the closed 
loop transfer function of the servo. 
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This development will parallel the Patton-Abbott approach so that 
from the same given data (closed loop transfer function, factored form) 
it is possible to obtain rapidly the transient response from their 
paper and the frequency response from this paper. This paper will 
discuss in sufficient detail the development of the frequency response 
approach so that no reference to the Patton-Abbott paper will be required. 



- 2 - 



1 . GENERAL DEVELOPMENT 



Since the time domain or transient response is used so widely in 
the initial phases of servo study, all developments in this paper 
originated with the closed loop transfer function expressed in Laplace 
notation (s-plane) with all initial conditions equal to zero. The 
transformed equation was then converted into frequency response notation. 
For this to be a valid conversion, certain basic definitions are in order, 
plus some additional restrictions and assumptions used in this devel- 
opment. The frequency response of a system was defined in part by 
considering only steady state conditions. Thus, all transient conditions 
have vanished. To convert the transfer function from Laplace notation 
into the frequency response equation, it was only necessary to make a 
substitution of the variable. This was based on the assumption that 
the input was a pure sinusoid and that the output was also a pure 
sinusoid since the steady state condition exists. Furthermore, this 
forces the restriction that the system must be linear. The complex 
variable s * o* -I- ju) then becomes the imaginary frequency variable 
since the real part equals zero (a = 0). Therefore s » ju); and when 
this variable substitution was made in the closed loop transfer 
function, the frequency response equation resulted. 

2 . 1 Basic Definitions and Systems Characteristics 

At this point, a tabulation of the various definitions, restric- 
tions, and assumptions would be in order: 

a. The system was linear. 

b. The system was stable. 

c. The system was characterized by unity feedback, i.e., 

M = 1 at 0) = 0. 

do The closed loop transfer function was known in its factored 

form. 

Based on the above, some general frequency response characteristics 
could be stated. 
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2 . 2 M- Con tour 



a. At zero frequency, M is unity. 

b. At zero frequency, the slope of the M contour is zero. 

c. As the frequency approaches infinity, M approaches zero. 

d. As the frequency approaches infinity, the slope of the M 
contour approaches zero. 

e. At the maximum or minimum value of M, the slope of the M 
contour is zero. 



2 . 3 N-Contour 

a. At zero frequency, N is zero. 

b. As the frequency approaches infinity, N approaches ^ tt 
degrees . 

(m is the number of zeros and n is the number of poles of the closed 
loop transfer function) 

2.4 Parameters 

To display graphically the frequency response, parameters were 
selected to define specific points on the M and N contours. The 
selected parameters were defined as: 

a. M : The value of M where the slope of the M contour was zero 

P 

(Does not include the value of M at zero frequency when M ■ 1.0). 

b. (Up* The value of omega (U) ■ 2Trf) where the slope o,f the M 
contour was zero. (Does not include cu ■ 0) . 



c. 


Np- 


The 


d. 


«bw' 


The 


e. 


“bw' 


The 


f. 


Nbw^ 


The 


g< 




The 


h. 


(U * 
1' 


The 


i. 


N^: 


The 



The bandwidth frequency (when M ■ .707). 



The notation when M ■ 1.0, except at cu ■ 0. 
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j. U)^: The undamped natural angular frequency. 

k. M : The value of M at the undamped natural angular frequency (00 ) 

oon n 

l. N • The value of N at the undamped natural angular frequency (00 ) 



The above selected parameters defined the general shape of the M 
and N contours. However, it was discovered that with only these points, 
the M and N contours were in some cases too general; the contours had 
insufficient shape to be of significant value. Better approximation of 
the M and N contours were provided by the selection of additional M, 
and tw points to *'f ill-in** the shape of the contours. It was discovered 
that, mathematically, the choice of as the basis for these additional 
points reduced the M and N equations to relatively simple relationships 
and insured that the additional points were located in the critical 
regions of the contours. The benefits gained by the additional points 
more than compensate for the slight additional labor by providing a much 
more accurate approximation to the M and N contours. The additional 
parameters selected were: 

m. An arbitrary frequency proportional to ((D^ ^ k ^p) > 

where k was a selected constant as discussed later. 



n . 




The 


va lue 


of 


M 


at 




o . 




The 


value 


of 


N 


at 





A sketch of the general M and N contours is shown in Fig. 2-1. 



2. 5 Classes of Servomechanisms Studied 

From Fig. 2-1, it can be seen that it was relatively easy to 
construct an approximate frequency response if the above defined para- 
meters were known. This paper contains the above parameters presented in 
graphical form. This paper also contains the technique utilized in the 
parameter determination such that the frequency response can be rapidly 
approximated for specific classes of systems. The three systems covered 
in this paper were: 

I. Pure Second Order - The closed loop transfer function consisted 
of two complex poles only. 
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Fig. 2-1 General M and N Contours Showing Parameters Chosen 
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II„ Second Order with One Zero ^ The closed loop transfer function 



consisted of one pair of complex poles and one real zero. 

III. Pure Third Order - The closed loop transfer function consisted 
of one pair of complex poles and one real pole, 

2 6 Presentation of Parameters 

The graphical presentations for the pure second order system were, 
for the most part, loci of constant values of the various parameters 
plotted on the s^plane. For the other two systems covered in this 
paper, it was necessary to express the coordinates of the complex poles 
as a function of the magnitude of the real pole or zero. In accomplishing 
this, the angular frequency, CJO, also became a function of the real pole 
or zero location. Thus the magnitude and phase shift for the frequency 
response was shown to be a function of the complex pole location relative 
to the real pole or zero. The graphical presentations for the second 
order plus one zero and for the pure third order systems were plotted in 
two ways. First, by loci of constant values of the various parameters 
plotted on a referred complex plane, that is, with ordinate and abscissa 
being the complex pole location relative to the real pole or zero. 

Second, by loci of constant referred imaginary part of the complex pole 
plotted against referred real part of the complex pole as abscissa versus 
the chosen parameter as ordinate, 

2 . 7 Digital Computer Utilization 

The solutions of the equations that will be presented later when 
a specific system is discussed were obtained from a Control Data Corp- 
oration high-speed digital computer Model 1604, located at the U. S. Naval 
Postgraduate School, Monterey, California (Appendix A). Without the 
availability of a high speed digital computer, it is doubtful that 
this paper could have been developed without the utilization of an 
analog computer plus excessive labor. 
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3 THE PURE SECOND ORDER SYSTEM 



The pure second order system treated in this section was pre- 
viously defined as a system whose closed loop transfer function consisted 
of only two complex conjugate poles. This system has been extensively 
treated in most servo textbooks, usually in terms of the damping ratio 
and the undamped natural frequency. Here, however, it will be treated 
by the rectangular location of the closed loop poles. The closed loop 
transfer function of the pure second order system was stated as /!/: 



f^Cs) 



'■) I'-i 



Cl 'f' b ^ 



(3.1) 



The poles of this system were plotted on the s-plane in Fig. 3-1. 
Converting Eq. (3.1) to the frequency response form, it became 




3 . 1 M and U) Calculations 

Equation (3.2) rearranged became 




from which 



(3.2) 



(3.3) 



M 






) 




I 






(3.4) 



Note that if the system gain were constant, then "a” and "b" would be 
fixed and the magnitude would be a function only of system frequency, 
i . e . , M = f ( CD) 
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Fig. 3*1 Poles of the Pure Second Order System (s-plane) 
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3,2 M and U) Determination 

J E 

CD was obtained by finding the value of CD where the slope of the 
P 

M versus CD curve was equal to zero« For ma thcna t leal convenience > this 
was done by first squaring Eq . (3.4), then taking the first derivative 
and setting the numerator equal to zero. 



(cl -» ,b ) (Z (a 4 - h ^Cu^)(~Zco) -f 8 o^co j — 0 



(3.5) 



Simplifying, this becomes 

•‘-CL 

and CD = 0 (see Section 2,4b) 



(3.6a) 

(3.6b) 



Substituting Eq . (3.6a) in Eq . 



(3.4) gave the value of M at CD 




Cl, D 

Zob 



(3.7) 



Note that from Eq. (3.6), M would not occur whenever **a” was less than 

P 

"b”, i.e., when the real part of the pole was greater than the imaginary 
part. (Figs. 3-2 and 3-3) 

3.3 CD, Determination 
bw 

The equation for the bandwidth frequency was obtained by substituting 
the bandwidth magnitude (M = .707) into Eq . (3.4) and solving for the 

frequency This resulted in 



(Fig. 3-4) 
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1 



3.4 M and ud Determination 
cJDn n 



The undamped natural frequency was defined as 



cOr, a. fh 

This equation expressed U)^ directly in terms of "a" and "b". 
The Substitution of into Eq. (3.4) yielded 



(3.9) 
(Fig. 3-5) 





(3.10) 



Equation ( 3. 10)expressed M in terms of the closed loop poles of the 

CJDn 

systemo (Fig 3-6) 

3 5 CJD^ Determination 

As defined previously, 00^ was the frequency, other than zero, at 
which M was unity. Letting M = loO in Eq . (3.4) and solving for 
yielded 



Lev, — 




(3.11) 



Simplifying this expression by recalling 

resulted in (JJ,^ ^2^ OU-p (3.12) 

(Fig. 3-7) 

3 6 m, and CJO, Determination 
k k 

As stated in Section 2, intermediate values of M and N, in addition 
to those determined above, would result in a smoother plot approximating 
the M and N contours. The values of k chosen for this class of servo 
were .5 and 1.5 These particular values were selected as being those 
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points which would be most beneficial in smoothing the plot of the 
M and N contours^ and because the equation was the same for k = .5 



as tor k = 1.5. 


In Section 2, 




CV^ CC k COp 


or 


\fk,ujp ^i~5cop 


and 


1.5 


(Fig 3*7) 




The value of 


and 0)^ were substituted into 










2 



! Z d '' h' 









(3.13) 



Where = Mj^ , recalling . Loci of constant values 

of were plotted in relation to "a" and "b". (Fig. 3-8) 



3 . 7 N Calculations 
From Eq. (3. 3) 



r A-. )_ - » A , „ 



t I 'X X 
. 4 0 -cO 



+Zai 



— ^ M 



O) 



from which, 



/V " -dtrciao 



0.1 



X , tX. X \ 

ti. 'th - iy J 



(3.14) 



Equation (3.14) was then the basic N equation for the pure second order 
system. 
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I 



3.8 N Determination 

_2 

0)^ was previously calculated by Eq . (3.6). The substitution of 
Eq o (3.6) into Eq« (3.14) resulted in the solution of related 
directly to "a" and "b”, the closed loop poles of the system. 




— o. ret an 




(3.15) 



Loci of constant values of N were then determined for all possible 

P 

"a*' and *‘b** values. (Fig. 3-9) 



It is pertinent to note that the equations presented were in most 
cases not solved as statedo Algebraic manipulations were performed to 
make the equation compatible to digital calculations by writing the 
equations in terms of "b*' as a function of *'a*‘ and the particular 
pa rameter . 



3.9 N Determination 
(Dn 



N was determined in a manner similar to N . The value of U) in 
(A)n p n 

terms of “a” and *'b" was stated in Eq . (3.10). This relationship of 
CJt)^ was substituted in Eq. (3.14) and resulted in 



AJ ~ circ'ta n 



Z 






L CA, 



4 



i- /- ri 
b -( = a 



■ht 



- 



which simplified to 




-arcian 



f Za\a-t h" -' ] 

L 0 1 



(3.16) 



or 



N = “90 
U)n 



Thus, for a pure second order system was a constant and equal to -90 
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3.10 N, Determination 
bw 

The equation for in terms of "a*' and "b" was obtained by 

substituting Eq . (3.8) into Eq . (3.14). This resulted in 



|\/, —-arctan 

dW 



L ZcL^- 1 .4:14 l/aZ-t ' -* 



(3.17) 



(Fig. 3-10) 



3.11 Determination 

Equation (3.12) determined the values of O) when M = 1.0 for non-zero 
frequencies. When the expression of CD^ , determined from Eq . (3.12), 
was substituted into Eq. (3.14), 






(Xf-Cl 



c'ta 



n 



Za, 









a 






- i-A 



(3.18) 



Solutions to Eq. (3.18) were plotted as constant loci for all values 
of *’a*' and *’b’' . (Fig. 3-11) 



3.12 N, Determination 
k 



As stated above, the k values selected were .5 and 1.5 where 

^ In a manner similar to the determination of above, 




(3.19) 



from which values were determined explicitly in terms of "a'* and *'b" 
for both values of k. (Figs. 3-12 and 3-13) 



3.13 Use of the Curves 

At this point, all the selected parameters necessary to sketch the 
M and N contours have been determined from the location of the closed 
loop poles. With the values of the closed loop poles known and the 
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curves plotted as described above, a rapid graphical approximation of 
the frequency response of the pure second order system can be constructed. 
For convenience, the following recommended procedure is presented for 
the use of the pure second order curves: 

a. Enter Fig. 3-2 with the complex pole coordinates -a + jb, and 

obtain M . Enter Fig. 3-3 with -a + jb, and obtain uo . (Note that U) 

P P P 

is the ordinate in Fig. 3-3.) 

b. Similarly obtain U), (M = .707) from Fig. 3-4, (Ji) and M 

bw n CA)n 

from Figs. 3-5 and 3-6 respectively. 

c. Obtain O)^ (M = 1.0) from Fig. 3-7, entering with the 
determined above. 



wi 



d. Obtain for both ^nd IC 2 = 1.5 from Fig. 3-8, entering 

th "a" and "b". Determine ^ ' 60^), 



and 



f W- - 



■\/r 



5 tOpJf 



from Fig. 3-7, entering with the value of 0)^ obtained from step (a) above. 

The M contour for the specified closed loop poles may then be 

sketched, having easily obtained the M and corresponding U) for six known 

points and realizing that M "= 0 and has a slope of zero at (A) = 0, and 

also that the slope is zero at U) and at U) = 

P 

The N contour points are obtained in the same manner by using 

Figs. 3-9 through 3-12 to obtain the values of N corresponding to the 

same frequencies found for the M contour above. The one exception is 

N , which is always equal to -90^ for this system. Using these N, 

Ci)n 

(JO values and knowing that N = 0 and has zero slope for (JO = 0 , and also 
that N = -180^ for Q) = for a pure second order system, the N contour 
can now be accurately sketched. 



Note that this paper does not present curves which will give any 

values at frequencies higher than oo , since it is unlikely that there 

bw 

would ever be any interest in the region beyond the bandwidth frequency. 
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3 1 4 Exarrple 

Given: *a ^ - -1 yl 

Solutdon: 

SPECIFIC DATA 

TABLE OF PARAMETER VALUES 



X 


0) 


CJL) Fig. No. 


M 


M Fig. No. 


N 


N Fig. No. 




X 


X 


X 


X 


X 


X 


p 


Ift 


3-2 


1.Z5' 


3-3 




3-9 


bw 




3-4 


. 707 


N.A. 




3-10 


0) 

n 


Z25" 


3-5 


1. 12 


3-6 


-90° 


N.A. 


1 




3-7 


1.0 


N.A. 


- 2 


3-11 


*^1 


172 


3-7 


11/6“' 


3-8 




3-12 


•^2 


2.10 


3-7 


l.i/^ 


3-8 


'82"^ 


3-13 



GENERAL DATA (for all pure second order systems) 



U) 


M 


dM/du) 


N 


dN/dcu 


0 


1.0 


0 


0° 


0 


U) 


M 


0 






p 


P 








00 


0 


0 


-180° 


0 



Fig 3-14 (page 17 ) compares the frequency response obtained from 
this rapid approximation method to the actual M and N contours obtained 
from a digital computer program of the general M and N equations. 

Note : Although a scale is shown on these pure second order curves for 

convenience, the loci are really independent of scale as they depend 
only on the ratio of *'a" to *'b” Hence, any convenient scale for O’, 
and may be substituted for the scales shown as long as all 
scales are in agreement 
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4. THE SECOND ORDER SYSTEM WITH ONE ZERO 



The addition of a real zero to a pure second order system 
resulted in a system whose closed loop transfer function was 




} 



d 



(s i-d) 



(4.1) 



The singularities of this system were plotted on the s-^plane in 

2 2 

Fig 4“1. The constant a + b , usually interpreted as system 

d 

gain, reflected the assumption of unity feedback as stated previously. 

By letting s = j(D, the frequency response form of Eq. (4.1) 
resulted in 






(4.2) 



This expanded and rearranged became 

db^)-4 j (ojg-huj 

(cla~i (zdaco) 





(4.3) 



Dividing the numerator and denominator of Eq . 



(4.3) by 




(4.4) 



From Eqo (4 4), it was seen tha t "aV'b'", and U) may all be referred to the 
zero location, d. 
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s- plane 








Fig. 4*1 Poles and Zero of the Second Order System with One Zero 
(s- plane) 
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From Eq . (4.4), 




(4.5) 




(4.6) 



and 



^ Ml 



(4.7) 



These quantities were illustrated in Fig. 4-2. Equation (4.4) then 
became 



Thus the frequency response for the second order system with one zero 
became a function of the "referred** location of the complex roots and 
the **referred** frequency. This means then that the frequency response 
was the same for all systems having the same relative pole-zero con- 
figuration when the frequency was also referred to the zero location. 
(Fig. 4-3) 

4. 1 M and Q Calculations 

Equation (4.8) was expressed in terms of its magnitude and phase 
angle by 
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m-plane 



-i 







X 












->• Wy 
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Figo 4-2 Referred Quantities of and w^, Second Order System 
with One Zero 
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System A: 



System B: 



a = 1.0 



a = .01 



b = 2.0 
d = 1.0 



b = .02 
d = .01 



"*1 ‘/d/ 



1.0 



m., 



/d/ 



1.0 



ra. 



■ /d/ 



m 



2 



b 

SI — 

/d/ 



2.0 




0 

0 



1 



.01 



2 3 

.02 .03 



4 (d * 1.0) 

.04 m (d = .01) 

B 



Fig. 4*3 Frequency Response of Two Systems with the Same Relative 

Pole*Zero Configuration, Second Order System with One Zero 
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From this, 







(4.10) 



Equation (4.10) was the basic M equation for the second order system 
with one zero. 



4.2 M and Q Determination 

-E. P 

The equation for was obtained by first squaring Eq. (4.10), 
then differentiating with respect to Q and setting the numerator of 
this equal to zero. This reduced to 






(4.11) 



Substituting this into Eq. (4.10) and simplifying, 



Mp= 






nrit -hX 7nf -h 






(4.12) 



Loci of constant values of m 2 were solved by digital computer and were 
plotted against versus M and Q . (Figs. 4-4 and 4-5). 

1 p p 



4.3 Q, Determination 

DW 

The equation for the referred bandwith frequency was obtained by 
substituting M ® .707 into Eq. (4^10) and solving for After 

considerable algebra, this reduced to 




(Fig. 4-6) 



Mf +yyi^ + 



{z(m^+yn^f [ 



7n, 




(4.13) 
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( 



4.4 M/-S and Q Determination 
L2n n 

The referred undamped natural frequency was defined for this 
system as 




1 






(4.14) 



This expressed Cl as a function of and m« . (Fig. 4-7) 
n 1 / 

Substituting Eq, (4.14) into Eq . (4.10) and simplifying yielded 





(Fig. 4-8) 



(4.15) 



4o5 Cl^ Determination 

As defined previously, Cl^ was the referred frequency other than 
zero at which M is unity. The equation for Cl^ was found by letting 
M = IcO in Eq. (4.10) and solving for Cl^. 










(4.16) 



(Fig. 4-9) 

4.6 and 

The use of the curves mentioned up to this point for this class 

of system provided the four main points required to sketch roughly the 

M contour. In order to '*fill-in** the blank portions of the contour on 

both sides of Q , it was decided to provide curves for three additional 
P 

points on the M contour in the same manner as was done for the pure 
second order system. The use of these additional points was not 
necessary to sketch the frequency response; however, their use provided 
a more accurate approximation of the M contour. 






Determina tion 
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As in the pure second order system, the following definition was made 






where k is a constant 



Substituting into Eq . (4 11) for , 



(4.17) 






(4.18) 



Values of .5, 1.5 and 2,5 were selected for k. M, solutions were 

k 

obtained for each k value by solving Eq< (4.18) for and then sub- 
stituting this value into the basic M equation with the same values of 
m^ and m^« This was accomplished handily by one digital program. 

(Figs 4-10, 4-11, 4-12,and 4-13) 

4. 7 N Calculations 

Referring to the basic closed loop transfer function in frequency 
form, Eq. (4.9), the general phase angle equation for the second order 
system with one zero, was written as 



W= arci 



Si (lyi ^ 4“ —XL ^ 



V yrS 



(4.19) 



The previously obtained parameter frequency equations (Q , Q, , ^ , 

p bw n 

, and the three ^’s) were solved by digital computer programs. The 
resulting parameter frequencies were then substituted into the general 
N equation (4.19) along with the corresponding values of m^ and m^ • The 
following figures resulted, all plotted with constant m^ loci against m^ 
versus N: 



Fig 4-14 N as function of ro- and m^ 

^ D 12 



Fig. 


4-15 




as 


function of 


m^ and 










Fig. 


4-16 


^Qn 


as 


function of 


ro^ and 


m2 








Fig. 


4-17 


as 


function of m^ and m^ 








Fig 


4- 18 


i<i 


as 


function of 


m^ and 


m^ 


for 




= .5 


Fig 


4-19 


k2 


as 


function of 


m^ and 




for 


^2 


= 1.5 


Fig. 


4-20 




as 


function of 


m^ and 


"'2 


for 


^^3 


= 2.5 
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4 , 8 Example 

Given: Complex poles; -a + jb = -2„0 + j 2.0 

Real zero; - d =“2.0 



Solution: 

SPECIFIC DATA 

TABLE OF PARAMETER VALUES 




GENERAL DATA (for all Second Order plus One Zero Systems) 



CJD 


M 


dM/du) 


N 


dN/duu 


0 


1.0 


0 


0° 




0) 


M 


0 




- 


P 


P 








00 


0 


0 


O 

o 


0 



Fig 4*21 (page 39 ) compares the frequency response obtained from 

this rapid approximation method to the actual M and N contours 
obtained from a digital computer solution of the general M and N 
equations o 
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5 THE PURE THIRD ORDER SYSTEM 



The addition of a real root to the pure second order system 
resulted in a new system which will be referred to as the pure third 
order system. The closed loop transfer function of the pure third 
order system had one pair of complex conjugate poles and one real pole. 

The closed loop transfer function for the unity feedback third 
order system was written as 



Letting s = juo in Eq (5.3) resulted in the frequency response form 
of the transfer function 




(5.1) 



Fig 5-1 defined these roots on the s-plane. 



Equation (5»1) was simplified as follows: 




(5.2) 



and expanded became 






iC f-'-i -f z tie i- 1')'} + c 



(5.3) 





^ ^ • Z c- f - 4 C ‘-y) (5.4) 



or 
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Fig. 5-1 Poles of the Pure Third Order System (s-plane) 
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Numerator and denominator of Eq . (5.5) divided by c : 




(5.6) 



n i< 

Observing that a power of c exists in all terms corresponding to 
the same power of **a** , **b** and (i); the following definitions were 
made : 





(5.7) 



^2 



A 




il 



A 40 




(5.8) 

(5.9) 



These quantities were illustrated in Fig. 5-2. Equation (5.6) then 
became 




7Cj 



2 - 






A 5^ 

= Me- 



(5.10) 



n^ and n^ were the complex pole locations referred to the real pole 
location and Cl was the referred frequency. The frequency response 
for the pure third order system was then a function of the referred 
complex pole and the referred frequency. Thus all systems having the 
same relative pole configuration will have identical response if the 
frequency were also referred to the real pole location. 



^ ^ M and Q Calculations 

The magnitude portion of Eq. (5.10), after algebraic manipulation, 
reduced to 




(5.11) 
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^ n^c *- 



-X 



s-plane 






n 



1 



/c/ 



n 



2 




X 



n-plane 




Fig, 5-2 Referred Quantities and n 2 , Pure Third Order System 
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Equation (5.11) was the basic M equation for the pure third order 
system. 



5 . 2 Discussion of the M Contours 

It was found that the M contour for a pure third order system 

will be one of three basic "shapes" as shown in Fig. 5-3a, b and c, 

depending on the values of n^ and n 2 . Figure 5- 3a was the situation 

when there was no peak. No value of Q was obtained and hence no values 

P 

for or the Figure 5-3b was similar to the second order M 

contour shape and had a single peak. For this general "shape" points 
were obtained for Sc Sc * 



The third general "shape" is Fig. 5* 3c, which resulted from the 

equation (5.12) which yielded two non-zero frequencies at which 

the M contour slope was zero, one being a "peak" frequency and the 

other being a "minimum" frequency. This latter frequency was designated 

n to avoid being confused with the Q . It was also discovered that 
m ° p 

with this shape of M contour, it was possible to have M less than 

P 

unity. Another possibility was to have less than and at the same 
time greater than one. This resulted in two values of two values 

of r^k > ^ ^ > and (one relative to and the other relative to 

M ), and three values of U (one for each crossing of the M = .707 line), 
m DW 

in addition to the usual points at Q , C , and Q . 

p m n 

Figure 5-4 might be thought of as the complex n-plane and shows the 
regions in which certain n^ and n 2 combinations produce one of the 
three above M contours. This figure indicates the general shape of the 
M contours that result from the known n^ and n 2 values. 



5.3 M and Q Determination 
2 

The equation for Q was obtained by squaring Eq . (5.11), then 
differentiating with respect to Q and setting the numerator of this 
equal to zero. This reduced to 
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Fig. 5-3 The Three Basic M Curve Shapes for Pure Third Order System 



(a) 



I 






r 



The plus and roirus values of the inner square root quantity were both 
valid and resulted in two possible solutions for The consequences 

of two solutions were discussed in detail in paragraph 5.2 above. 

Loci of constant values of n^ were solved by a digital computer 

program and were plotted against n_ versus Q and Q in Figs. 5-5a 

1 pm 

and 5“5b In the same computer program, values of were obtained 

for corresponding values of n^, and n^ by substituting these 

known quantities into the basic M equation, Eq . (5.11). Figures 5-6a, 

5^6b, and 5~6c contain M plotted as a function of n_ and n^ . No 

P 1 ^ 

specific M equation was derived due to the obvious complexity 
involved in substituting Eq. (5.12) into Eq. (5.11) for Q, 



5.4 fl Determination 

DW 

The equation for the referred bandwidth frequency was obtained 
by substituting M = .707 into Eq (5.11) and solving for n^ as a function 

2 



of Cl and n 

DW 



This involved extensive algebra which reduced to 






■al, -tKJ-zrCn)-, 






(5.13) 



The + signs were both valid here also and their effect will be discussed 
later. (Fig. 5-7) 

5.5 and Q Determination 

L2n n 

The undamped natural frequency for the pure third order system was 
defined as 






“-3 



-a “ 

Cl^^ was expressed as a function of n^ and n 2 



(5.14) 



(Fig. 5-8) 



Substituting this into Eq. (5.11) and simplifying yielded 










(5.15) 



MCxi thus expressed as a function of n^ and n 2 * (Figs. 5-9a and 5-9b) 

5.6 Determination 

The referred frequency, other than zero, at which M “ 1.0 was 
defined as Q^. 

Letting M * 1.0 in Eq. (5.11) and solving for resulted in 



i)., 






/l,4 



^ I- 

- .5 



±V^ 



7^1 -"•471( 71^ 'i~. 



X 



Z5 



(5.16) 



(Figs. 5-10 a, b, c, and d) 

5.7 and Determination 

For the same reasons as discussed in Section 2.4, it was found 
desirable to obtain "fill-in" points on both sides of 0^, This 
was done by defining 

- kOp (5.17) 

which eliminated the square root of k as used for the previous two 
classes of systems. 

The values of k chosen were .8, 1.4, and 1.8; therefore 

fl * . 8 fip 



(1 ■ 1 . 4 Qp 

k2 



Q, “ 1 . 8 Cp 

k-3 



(5.18) 

(5.19) 

(5.20) 
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Fig. 5-11 was the almost trivial plot of these three equations, 
included for optional use in determining the values once has 

been determined. 

Loci for the three were calculated by a digital program which 

utilized Eqs. (5.11), (5.12), (5.18), (5.19), and (5.20); thus it was 
not necessary to attempt to obtain any specific equations. 

(Figs. 5-12 a,b, 5-13 a,b, 5-14) 



5 . 8 Determination of N Curves 

Referring to the basic closed loop transfer function in frequency 
form, Eq . (5.10), the general phase angle equation for the pure third 
order system was written as 



aKtan f +2.n,+nf,) 



(5.21) 



The previously obtained parameter frequency equations (^^j 
, and n^) were solved by digital computer programs. The resulting 
parameter frequencies were then substituted into the general N equation, 
Eq. (5c 21), along with the corresponding values of n^ and 1 X 2 * The 
following figures resulted, all plotted with constant n^ loci against 
n^ versus N: 



Fig. 5-15 a 
Fig. 5-15 b 
Fig. 5-16 a 
Fig. 5-16 b 
Fig 5-17 a 
Fig. 5-17 b 
Fig. 5-17 c 
Fig. 5-18 a 



N as a function of n^ and n^ 

p 1 2 

N as a function of n^ and n^ 

m 12 

N, as a function of n^ and n^ 
bw 12 

Nbw ^ function of n^ and n 

as a function of n^ and r \2 
(expanded) as a function of n^ and n 2 
(2nd value)as a function of n^^ and n 2 
for maximum (^p) 
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Fig. 


5-18 


b 




for 


minimum 




Fig. 


5-19 


a 


Nk 

k2 


for 


maximum 


<%) 


Fig. 


5-19 


b 


‘^2 


for 


minimum 




Fig. 


5-20 






for 


maximum 


<%) 



An equation for was obtained by substituting Eq. (5.14) into 
the basic N equation, Eq . (5.21), 



— o^rciart j -j 



(5.22) 



Loci of were plotted as a function of n^^ and n 2 • (Fig. 5-21) 
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5 . 9 Example 



Given Complex poles* i J 



Real pole: 



-c =-2.0 



Solution : 
SPECIFIC DATA 



"l = /c/ = 



TABLE OF PARAMETER VALUES 



X 


Cl 

X 


(JU k'k'k 

X 


C. Fig. No. 

X 


M 

X 


M Fig. No. 

X 


N 

X 


N Fig. No. 

X 


p 


/ 8S 


■3-20 


5- 5a 


//^ 


5-6a ,b 


~J3C" 


5 - 15a 


m 




/.i3^ 


5- 5b 




5-6c 




5-15b 


bw 


3^ 


6 c 


5-7 


. 707 


N.A. 




5-16b 


bw 


N A. 


— — 


— 


. 707 


N.A. 


N.A. 


— 


bw 
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GENERAL DATA (for all pure third order) 
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Fig. 5-22 (page 69 ) compares the frequency response obtained from 
this rapid approximation method to the actual M and N contours 
obtained from a digital computer program of the general M and N 
equations . 
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6. CONCLUSIONS AND RECOMMENDATIONS 



The following conclusions were made based upon the investigation 
of the three classes of servomechanisms studied in this thesis: 

a. All pure second order servomechanisms with the same relative 
complex pole location will have the same frequency response. 

b. For the pure second order servomechanism, the parameter loci 
plotted are made applicable to any location of complex poles in the 
left half of the s-^plane by changing the scales on the complex plane 
and on the parameter frequency loci. 

c. The frequency response for servomechanisms with three 
singularities can be expressed as a function of the referred real and 
imaginary parts of the complex poles and the referred frequency where: 

1 The referred real part of the complex poles is the 
actual real part of the complex poles divided by the 
magnitude of the real pole or zero. 

2. The referred imaginary part of the complex poles is the 
positive imaginary part of the actual complex pole 
divided by the magnitude of the real pole or zero. 

3. The referred frequency is the actual frequency divided 
by the magnitude of the real pole or zero. 

d. All servomechanisms with three singularities which have the 
same re lative pole -zero locations will have the same frequency response 
relative to the referred frequency. 

e. Curves for determining values of several optional *’fill-in** 

parameters were included in addition to the normally accepted parameters 

of M ,0) and O)^ . 
pO) p bw 

f. An accurate graphical approximation of the frequency response 
of a servomechanism can be rapidly obtained using the curves produced 
in this thesis- 



It is recommended that the following areas be considered for further 
investigation using the technique presented in this paper; 

a- Servomechanisms with more than three singularities, and 
b. Systems other than servomechanisms. 
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8 . Appendix A 



The Control Data Corporation Model 1604 digital computer 
important characteristics are summarized as follows: 

1. All transistorized 

2. 48 bit word length 

3 Stored program 

4 Single-address logic, two instructions per 48-bit word 
5. Indirect addressing feature 

6 32,768 48-bit words of magnetic core storage: 

(a) Stored in two independent 16,384 word banks, 
alternate phased 

(b) 48 micro ' seconds effective cycle time 
(representative program). 

(c) 6.4 microseconds total cycle time. 
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